How Sensitive is the CMB to a Single Lens? 



Ben Rathau£], Anastasia FialkoA|l| and Nissan Itzhakj§ 

Raymond and Beverly Sackler Faculty of Exact Sciences, 

School of Physics and Astronomy, Tel- Aviv University, 

Q ' Ramat-Aviv, 69978, Israel 

(N 

! Abstract 

(N 

We study the imprints of a single lens, that breaks statistical isotropy, on the 

O I CMB and calculate the signal to noise ratio (S/N) for its detection. We emphasize 
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1 Introduction 



The role weak lensing plays in CMBphysics is well understood in models that assume 



statistical isotropy, such as ACDM [ifQ. Statistical isotropy assures that no preferred 
direction exists, and weak lensing acts to smear the peaks and troughs of the power 
spectrum by convolving different scales. 

There are, however, some unexpected features in the current data that, at least naively, 
appear to hint that statistical isotropy is broken at the 2-3cr confidence level. These include 
the drop in the power spectrum at large scales that is highly sensitive to the galactic 
plane mask 0-3], the alignment of the low / modes (sUl, the anomalously large bulk 
flow 12- 



14], the cold spot 15 -19] and giant rings [2^ (for recent reviews see [sl, [2l[ 22|). 
Cross-correlating these features with weak lensing of the CMB could provide us with some 
valuable information about the origin of these features. In particular, it was argued that 
a possible cosmological explanation to some of these features could be an anomalously 
large structure (e.g. 



18 



23 



26|). Such a structure, irrespective of its origin (which 



could be new physics or a statistical fluke within ACDM), may also leave a statistically 
significant imprint on the CMB via weak lensing. 

With this motivation in mind we study the CMB weak lensing signal associated with 
a single lens that breaks statistical isotropy. Our main goal is to determine under what 
circumstances will such a lens leave a statistically significant imprint on the CMB, that 



27 



is, a S/N larger than 1 (for previous works on the subject see 

The paper is organized as follows. In section 2 we summarize some known results 
about S/N calculations in two cases that are relevant for weak lensing of the CMB, 
assuming Gaussian distribution. In section 3 we first show that there is an upper limit 
on the S/N associated with a single lens. We refer to this upper limit as the ideal S/N. 
We also calculate the reahstic S/N from the temperature CMB data. Section 4 begins 
with a puzzle: we show that the expressions we got for the realistic S/N imply that 
an experiment with a resolution /^nax > 2000 has a larger S/N than the ideal S/N. We 
show that non-Gaussianities induced by ACDM weak lensing resolve this puzzle in an 
interesting way. Finally, in section 5 we calculate the S/N for two particular examples of 
large scale structure: a texture and a void. Both are possible cosmological explanations 
to the cold spot. We find that a texture cannot be detected via weak lensing, whereas a 
void is barely detectable. 



2 



2 Preliminaries 



In this section we review two deformations of a Gaussian distribution and calculate the 
associated S/N. As we shall see in the next sections both deformations are relevant for 
weak lensing of the CMB. 

The starting point is a system of n Gaussian fields with vanishing mean values 



(x*) = 0, i = l..n 

and non-trivial covariance matrix 

(xV) = C^^ 



The likelihood function associated with this system is 30N32| 



For simplicity we take Cq to be a diagonal matrix. 

Generally speaking, there are two classes of deformations of this system that do not in- 
duce non-Gaussianity. Since these deformations are from one Gaussian system to another 
Gaussian system the S/N associated with them can be calculated exactly. 

Class 1: Deformation of the mean value 

The Gaussian fields are shifted by constants (i.e. non-random variables) 

= x' + U. (2.2) 
This deformation modifies the mean values while leaving the covariance matrix intact 

{x^) = h\ {{x^ - {ms^, - {x^})) = Cl (2.3) 

The likelihood function of the deformed system is 

£(b) = , \ exp f--(x + b)TCn'(x + b)Y (2.4) 

^ ' (27r)"/2yd^ 2^ ^ ° ^ 7 ^ ^ 

To calculate the S/N associated with this deformation we proceed in the familiar way, 



using Fisher information theory (e.g. [30|-|32|). We define 

5(b) = -2 [log(£(b)) - log(£o)] = (x^Co 'b + b^Co 'x + b^Co 'b) , (2.5) 
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which is nothing but the |30|, |32|, and calculate its mean value which is the desired 
S/N ratio 

I) ' (b) = (5(b)) = I dx5(b)£o = b^Co ^b. (2.6) 

This class of deformations is relevant when searching for a known template, b, in the 
data. 

Class 2: Deformation of the covariance matrix 

We leave the mean values intact, (x*) = 0, and deform the covariance matrix 

Co ^ C. (2.7) 

In general, the deformed covariance matrix C is non-diagonal. The relevant likelihood 
function reads 

1 ^ ^ .Tn-U 



C(C) = ; exp — x^C-^x . (2.8) 

^ ' (27r)«/2yd^ 2 ; ^ ' 

To calculate the S/N we proceed in the same fashion as before. We define 

5{C) = -2 [log(£(C)) - \og{Co)] = x^ (C-i - Co ^) x + log (det C/ det Co) , (2.9) 
and calculate its mean value 

= (5(C)) = y"dx(5(C)/:o = Tr(CoC-^-l) +log(detC/detCo). (2.10) 

As is often the case, weak lensing included, one can expand C in some small expansion 
parameter e as 



C = Co + eCi + -C2 + ... (2.11) 

Expanding (12.101) in e we find that the linear term vanishes due to cancellation between 
the two terms in the r.h.s. of f l2.10p . This follows from the fact that C = Co is a 
local minimum. What is somewhat surprising is that the leading term, that scales as e^, 
depends only on Ci (and not on C2) 

Q \ 2 2 I Cfij 1 2 

^1 (2.12) 
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A familiar situation, which is irrelevant to our case, is that of diagonal Ci. Then the 
leading S/N reads 

1 V- f5G ^ ^ 



N / 2 ^ V c: 



with 5Cii = eCi. 

Applying this result to the Gaussian statistically isotropic CMB temperature (where 
i runs over / and m) we get the familiar expression 



where AC? = Ci^2/{21 + 1] 



3 Ideal and Realistic Expressions for S/N 

In this section we apply the deformations discussed above to calculate the S/N associated 
with weak lensing of the CMB induced by a single lens. We perform two calculations. One 
is of an ideal S/N in which the signal of the weak lensing by the single lens is competing 
only with the noise associated with ACDM weak lensing. The second is a realistic S/N 
in which the signal is extracted from the temperature data. In such a case, the signal 
is competing with various contributions to the noise on top of the noise associated with 
ACDM weak lensing. Thus we expect the ideal S/N to serve as an upper bound on the 
realistic S/N. 

3.1 Ideal S/N 

In ACDM the deflection potential field ip^{h) is statistically isotropic and follows a Gaus- 
sian distribution l|-3| 

(^L) = 0, {^^L^L') = Su'S^m'Cf, (3.1) 

where as usual we expand ip^{n) = Xl/m '^'m^m(n) and the angle brackets stand for an 
ensemble average. 

The effect of a single lens on ip(^ is to shift it by a constant Sipim that depends on the 
details of the single lens 

i^L^i^L + ^i^im. (3.2) 
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Since the shift is by a constant (and not a random field) the deformation associated with 
adding a single lens to ACDM is of the first class of deformations discussed in the previous 
section, and the ideal S/N is 



/ Ideal 



In the rest of the paper we work mostly in the fiat-sky approximation. The full- and 
fiat-sky variables are related in the following way |l| 



C{1) = Cu m = V^E^""^^^-^^'"''- (3.4) 

m 

In fiat-sky variables (13. ip takes the form 

(^^(1)> = 0, (^^*(1)^^1')> = (27r)2<5(l - V)C^i\), (3.5) 
and the ideal S/N is 



SV f d\ \6ij{\)f 



(3.6) 



where 1 is a 2D momentum. 

3.2 Realistic S/N 

The effect of the single lens on the CMB temperature is to shift the temperature at each 
point in a way that is determined by the defiection potential it induces, Sip- Expanding 
to leading orders we have 

TSL(n) = T(n) + V'6^ViT{h) + ^V'6ijV^^ViVjT{h) . . . (3.7) 

where T is the ACDM— lensed temperature field and T^^ is further lensed by the single 
lens. 

To calculate the S/N it is natural to define 6T{n) = T^^{n) — T(n) and treat it as a 
deformation of the first class. In that case we can use (12. 6p . to find 

I)' ={ST^C-^6T). (3.8) 

/ Temp 
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This is the point of view taken in [27| which in the flat-sky approximation yields for a 
spherically symmetric single len^ 



where 



/HI' 
^|r.(i-r)Wrc(|i-r|). (s.io) 

This seems correct since a known structure determines Sip and so it determines ST{h). 
Hence it appears that once again we are searching for a known template, this time ST{n), 
in the data. 

However, as emphasized in the previous section this approach can be used only when 
the template is a known constant deformation and not a random field. In our case 6ip is a 
constant deformation but (5T(n) is a random field. In particular, its mean value vanishes 
(5T(n)) = 0. We conclude that from the point of view of temperature fluctuations weak 
lensing of the CMB by a single lens is not a deformation of the first class. Hence (12. 6p 
is not a good starting point for the calculation of the S/N. In fact, it is a deformation of 
the second class. 

To see this and calculate the S/N it is convenient to Fourier transform (13. 7p and get 
the familiar result 

tsl(i) = T(i)-y ^r.(i-io#(r)T(i-n 

/ (l!y^ (1^^' ■ + ^" " ^] ^' ■ ^"nmin^nr + r - 1). (s.n) 

It is easy to see that indeed (T^'"(l)) = and that the covariance matrix 

Cov(li,l2) = (T^^*(li)T^^(l2)>, (3.12) 

is deformed. Eq. (I2.12p implies that to calculate the S/N to second order in Sip we need 
to know Cov(li, I2) only to first order. This is easy to calculate. For convenience we first 
define 

7(li, I2) = (li + I2) ■ IhCih) + 12^(12)] (3.13) 



^With the substitution / J2i 
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so that the off-diagonal terms of the covariance matrix are 

Cov(li,l2) = 7(-li,l2)5^(l2-li). (3.14) 

The diagonal terms of the deformation vanisl]0 so that the diagonal terms of the complete 
covariance matrix are just the original C(/)'s. 
Thus using fl2.12p we find 

SV _ 1 /• dl dr \j{-\,Y)5i;{V-\)f 
Ny'xen.p 2 7 (27r)M27r)2 C{l)Cil') " ^ " ^ 

In the case of a spherically symmetric deflection potential, for which = one 

can write fl3.15p as 



- 



where 

^'"/ (^^(iyf [i'-(i^(0 + (i'-iM|i'-i|))f • (3.17) 

Note that while Eq. fl3.17p looks similar to Eq. fl3.10p . because of the cancelation between 
the two terms in f l3.17p it is in fact smaller by roughly a factor of 10^. This makes a single 
lens much harder to detect via weak lensing of the CMB than expected. 



4 The Role of Non-Gaussianities 

In the previous section we have calculated both the ideal and realistic S/N associated 
with a single lens. As we now show it turns out that the realistic S/N is larger than the 
ideal S/N for experiments, such as Planck, SPT and ACT, with 6' resolution or better. 
This puzzle is illustrated in the context of a toy example in the next subsection. In the 
second subsection we show how non-Gaussianities induced by weak lensing of the CMB 
resolve this puzzle. 

4.1 The Toy-Model and the Puzzle 

So far we have had two types of signal to noise ratios: the ideal f l3.6p . and the realistic 
f l3.15p . We have argued that the ideal S/N must be an upper bound on any realistic S/N, 

^This statement is correct only in the flat-sky approximation. In full-sky, there is a small diagonal 
contribution to the covariance matrix, as we show in appendix [A] 
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and in particular on fl3.15p . This should hold for any single lens. Consider for simplicity 
a toy-model in which all the contributions to the deflection potential comes from a single 
mode lo 

<5V'(l) = (27r)2^5(l-lo). (4.1) 



The ideal S/N is then simply 



and the reahstic S/N is 

2 A ! I /■^ \ p J^ 1/11 im2 



(4.2) 



fsy ^ ASijjio) r dl i7(i,io-i)r 

VN/Temp 2 J {2nyC{l)C{\lo-\ 



Temp 

2 J (27r)2C(/)C(|lo-l|)' 
In light of the above argument we expect the ratio 



(4.3) 



p., . _ (S/N)Ten.p _ g-^(M / dl |7(M0-l)r 

2 y (2vr)2C(/)C(|lo-l|) 



to be smaller than 1 for all Zq- 
Writing the integral in f l4.3p as 

dl 



^ tS. ^d/ (4.5) 



(27r)2 7 (27r) 

quantifies the resolution capabilities of a given experiment via the cutoff LmaxE Calculat- 
ing S/Nxemp we find that for any Iq it scales quadratically with Lmax- Thus -R(/o) crosses 
1 for a sufficiently large Lmax, which is typically of order 2000 (see Fig. [T]). The quadratic 
dependence on L^ax follows from the fact that at large / the dependence of (S/N)Temp 
on C{1) drops and we are left with an integral over a constant. Therefore, this puzzle is 
generic and does not depend on the details of the model and/or parameters of ACDM. 

4.2 Non-Gaussianities Resolve this Issue 

As is often the case with puzzles, the origin of this puzzle is hidden in the assumptions 
we have made. In deriving (12.121) we assumed that the background temperature field 
follows a Gaussian distribution. However, it is well known that weak lensing induces 

■^Alternatively one can add the instrumental noise Ni ^ exp(P / 1'^^^) . 
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(a) 



(b) 



Figure 1: S/Nxemp for the toy model dijj oc 6{\ — Iq), as described in (14. ip through fl4.4p 
for a few values of Iq, normalized according to the S/Nideai (red dashed line). Here (a) 
shows the cumulative S/N as a function of Lmax for a wide range of Lmax while (b) 
a zoom- in picture of the cyan rectangle, emphasizing the region of L^^x where S/Nxemp 

crosses S/Nideal- 



IS 



non-Gaussianity [2|, |33[. In ACDM this non-Gaussianity is expected to start affecting the 



Gaussian behavior at / ~ 1000. At / ~ 2000 it modifies 



;he Gaussian results by about 5% 



33|. In fact, very recently this effect was measured |34l. |35| . and an agrement with theory 
was found. As we show below, when statistical isotropy is broken, the non-Gaussianity 
induced by weak lensing is more significant: It alters our findings by a factor of order 
unity exactly where it should, / ~ 2000, in order to solve the puzzle we have encountered. 

It is constructive to address this calculation via Feynmann diagrams. Let us first 
repeat the S/N calculation in the Gaussian case of the previous section using Feynmann 
diagrams. The Feymann rules associated with this Gaussian theory are as follows: 

• Assign a 2D momentum to each leg in the diagram, and a propagator, C{1), which 
corresponds to the ACDM (lensed) temperature power spectrum. 

• The single lens induces a two-leg "interaction" vertex 7(li,l2)/2, where 7(11,12) = 

7(ll,l2)5^(ll + l2)/C(/l)C(/2) 

^ 7(li,l2)/2. 



This vertex mixes between two momentum modes. 
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• Multiply each diagram by a proper symmetry factor. The relevant symmetry factors, 
and their derivation are elaborated on in Appendix [Bl 

The Gaussian calculation of the previous section is presented by the diagram of Fig. [2al 
Namely, it is a vacuum energy diagram with an external background field (that breaks 
statistical isotropy) induced by the single lens. Following the Feynmann rules we get, 
using ( I3.13P and the fact that C{1) and ip{n) are real, 

which is identical to f l3.15p . 

ACDM weak lensing induces non-Gaussianities in the temperature field. Assuming 
negligible correlations between the deflection potential and the primordial temperature the 
leading non-Gaussianity is a four-leg vertex. Hence we have to supplement our Feynmann 
rules with 

• ACDM weak lensing induces a four-leg vertex 

h I, 

= f(li,l2,l3,l4)/4! 



Since the contributions to the S/N we are interested in are quadratic in Sip, the relevant 



leading Feynmann diagrams are those of Fig|2bl In standard field theories (e.g. (p^) such 
diagrams are easy to calculate. Here, however, the vertex is rather complicated. The 
reason is that it has substructure that follows from the relation between the unlensed 
and the lensed correlation functions in ACDM, leading to topologically distinct diagrams. 
The details of the calculation are left to Appendix |Bl Here we merely quote the results. 
There are four different diagrams which contribute to the 2-loop S/Nxemp- One of these is 
positive, and the other three are negative. The sign of the overall leading non-Gaussian 
contribution to S/Nxemp is negative. Therefore, it lowers (S/N)|gj^p, as is illustrated in 

Fig. m 

Fig. [3]also shows that the leading non-Gaussian contribution becomes significant at / ~ 
900 and that for 1000 < /max < 1500 there is an approximated plateau in the accumulated 



^Leading here refers to the number of loops in a diagram or, alternatively, the nmriber of times the 
4-leg vertex appears in a diagram. 
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(a) 



(b) 



Figure 2: Visualizing the S/Nxemp as vacuum-energy Feynmann diagrams up to 2-loop 
order following the Feynman rules of section HI The diagrams here are schematic, and do 
not account for some complication due to ACDM four-leg vertex substructure, (a) The 
one-loop (Gaussian) part of the S/N. (b) The leading non-Gaussian contribution. 



S/N. The value of (S/N)^ at this plateau is about 1/10 of the ideal (S/N)^. These features 
are not sensitive to Iq. For /max > 1500 the accumulated S/N starts to drop. This is a 
nonphysical artifact of the fact that we keep only the first non-Gaussian correction. We 
expect higher (in loop counting sense) non-Gaussian corrections to remove this artifact 
and to alter the S/N for /max > 1500. It is beyond the scope of this paper to calculate these 
corrections and to determine the exact dependence of S/Nxemp on /max- We emphasize 
that, irrespective of /max, one cannot exceed ~ -^S/Nideai using the temperature data 
alone without performing this higher loop calculation. 

Physical situations, like the ones discussed in the next section, can be viewed as a 
superposition of toy models with different /q's. Therefore we expect them to behave in a 
similar way, namely 

WTen^p^TT^U 

to be a good approximation. 

5 Physical Examples 

In this section we study the weak lensing S /N associated with two examples of spherically 
symmetric structures that were proposed as possible cosmological explanations to the 
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(a) 



(b) 



Figure 3: The significance of tlie 2-loop contribution to tlie S/Nxemp- Here (a) sliows 
tlie 1-loop (blue) vs. (1 + 2)— loop (cyan) contribution to the integrand per /, while [(b) 
shows the (S/N)j^gj^j (red dashed), Gaussian part of (S/N)^gj^p (blue solid), and (S/N)jjg^[ 
(cyan solid and dashed) to leading non-Gaussian contribution for Iq = 50. The cyan line 
becomes dashed at the point where the overall non-Gassian contribution to the integrand 
becomes negative. 



WMAP cold spot (isl y, 1271, l28|. One is a cosmic texture (isi |23|, |27| and the other is a 



localized void (similar to that discussed in 



24 



27|,l28|). 



5.1 Cosmic Texture 



The deflection profile of a texture as discussed in [27|], is 

a{e) = At ^ 



(5.1) 



which is valid in the range < < Rt- We follow [27], and take Rt = 5° and At = 
5.19 X 10~^ to best describe the WMAP cold spot. In the flat-sky approximation we have 



27| 



TT 



5^(/) = -^^ArRre-^^^^^lRT + 2)//^ 
from which we can calculate the ideal S/N to find (using (13. 6p ) 



N 



/sky 2 X 10- 



(5.2) 



(5.3) 



Ideal 



where fsky is the observed fraction of the sky. We conclude that the texture's imprint on 
the CMB via weak lensing is too small to be detected. 
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5.2 Void 

The void case is somewhat more elaborate because there is a subtlety in the ideal S/N 
calculation that should be dealt with. Our starting point is the void discussed in |27| 
which was taken to have a sharp wall at Ry = 3.1°. Namely, the energy density of the 
void is 

(5.4) 



5p{d) [0 e<Rv 



Po I 1 9 > Ry ■ 

The deflection angle associated with such a void is 

\ Ave e<Rv 

with Ay = 0.01785. We see that the discontinuity in 6p induces a discontinuity in the 
flrst derivative of av{0) at 9 = Ry. 

Upon Fourier transforming (15. 5p . this discontinuity in position space translates to a 



signal in momentum space which does not decay 
Indeed, at large /, the deflection potential as in 



ast enough for the ideal S/N to converge. 



27| 



5i) = 47iAvRvMlRv)/l^ (5.6) 

scales like ~ /^^/^ while the deflection potential power spectrum of ACDM drops faster 
than so the S/N of the void diverges as anticipated (see Fig. Hj). This is a rather mild 
divergence (only logarithmic in /) that is not going to affect any calculation relevant for 
actual experiments. Still, this issue is worth addressing. 

The origin of this divergence is the assumption that the void is spherically symmetric 
with a sharp wall. In physical situations it is often the case that the thickness of the wall 
is much smaller than the size of the void and so a sharp wall is a good approximation. 
However, we do not expect voids to be approximately spherically symmetric. Hence 
a discontinuity in the flrst derivative of av{0) is very unlikely to happen. We have 
considered various smoothings of this discontinuity. Since the original divergence is only 
logarithmic in I we do not expect the details of the smoothing to considerably affect the 
S/N. Indeed, this is what we have found with a typical value of 



N 



~ 15 /sky. (5.7) 

Ideal 
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Sharp-wall void 



thed 



.Id 











100 



200 



300 



400 



'max 



Figure 4: The ideal S/N of the void, in two "versions": a void which is modelled by a 
deflection profile with a sharp wall (blue) as in (15. 5 p and which is seen to moderately 
diverge with Lmax, and a void whose profile was smoothed so that the first derivative of 
the deflection profile is continuous (red), which converges. 



Hence using (14. 7p the realistic S/N is approximately 



We conclude that in experiments with partial sky coverage and high angular resolution, 
such as ACT and SPT, /sky renders the void undetectable. Moreover, even in full-sky 
high resolution experiments such a void is barely detectable. 
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A Full-sky S /N 



The full-sky expression for the S/N in the case of a spherically symmetric lens reads 

2 

' ''"'^ (A.1) 



= E 

Temp /<;2,m 



Zm;Z2m 



ClCu 



2 '''^^ 



where C^^.i^^^ is the leading order correction to the covariance matrix, that is, 



Im;l2m2 



T* n^^^ 



^Im ^ hm2 / 1 



(A.2) 



and where Tim and are the /, m— harmonic moments of the ACDM lensed temperature 
field, and its first order correction due to the single lens respectively. 

The expression (lA.ip includes both diagonal and off-diagonal contributions, and im- 
plies that the covariance matrix is diagonal in the "m— dimension" (since we consider a 
spherically symmetric lens). The off-diagonal terms Ci^-um (with I ^ h), are derived in 



291] (Eq. 13), and here we mainly examine the diagonal contribution to the S/N (lA.ll) . 



A.l Diagonal Terms of the Covariance Matrix 



The diagonal terms are obtained by taking I = I2 in (IA.2P (or in Eq. (13) of 29|]). This 
gives 



(1) 

Imdm 



ciY^{-iri"{i" + 1) 



X 



/ / I" 





-m m 



121" + 1 
X (2/ + l)-^/^— 5V^po, 



(A.3) 



where 



m m' m" 

is the Wigner-3j symbol, and where we have fixed the coordinate system so that only 
m = modes of the deflection potential are non-vanishing. Hence, 



.(1) 

Imdm 



E 



(2/ 1) 
An 



I I I" 




[l"{l" + l)fC^5ijl,. (A.4) 
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We see that, unlike in the flat-sky approximation, the full-sky expression contains a non- 
vanishing contribution. 

Despite the fact that there is a diagonal contribution to S /N at first order we show now 
that the first order contribution to the estimator of the lensed power-spectrum vanishes. 
To see this note that since 



G 



lm,;lm, 

21 + 1' 

rri 

where X denotes the estimator of X, and since 



—m m 



(A.5) 



oc Si" 



(A.6) 



we get 



as claimed. 



(1) 



0, 



(A.7) 



A.2 S/N 



Collecting 
culated in 



3oth contributions, the diagonal of the previous section and off-diagonal cal- 
we find that 



where 



^/ Temp 



+ 



DP 



N 



ODP 



DP 



E 

Im 



Irri 



2Cf 



and where 



ODP 



4 ^ 

l<l2,l" 



I k I" 





(A. 



(A.9) 



(2/ + l)(2/2 + l) 

X {\l{l + l)-h{l2 + l) + l"{l" + l)]Ci 

+ Mk + 1) - l{l + 1) + /"(/" + 1)] Ci,f 6^f„o. (A.IO) 
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B Performing 2-Loop Calculations 



In section H] we have already outlined the Feynmann rules, which make the 2-loop calcu- 
lation easier. In this appendix we write down explicitly the "vertex-functions" 7 and F, 
draw all the diagrams that contribute to the S/N up to 2-loop order, find their symmetry 
factors, and then use the Feynmann rules to evaluate the different contributions. We then 
apply the toy-model deflection (14. ip to the results, thus simplifying them further. These 
results give together the 2-loop negative contribution seen in Fig. |3l All the calculations 
in this appendix are carried out in the flat-sky approximation. 

B.l The Vertex- Functions 

A vertex-function is the correlation function (in Fourier-space) divided by the incoming 
propagators. Therefore the 2-leg vertex function of a single lens is 

^^'''^ " Cih)Cil2) 

C{h)C{l2) 

The 4-leg vertex induced by the ACDM-weak-lensing effects is more complicated since it 
has non-trivial substructure. This follows from the fact that the vertex interaction of the 
lensed temperature fleld is calculated from the interaction of the unlensed temperature 
fleld and ip^, namely, 

'x; - ><: 

ii h ii 

where, all the legs on the l.h.s. are of the lensed temperature, and on the r.h.s. they are 
of the unlensed temperatures. The dashed line is the ijj^ propagator, C^. 
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For the above drawn orientation of the momenta the correlation function may be 
written as 



orientation 

(27r)' 6{h + I2 + I3 + I4) I3 • (li + I3) I4 ■ (I2 + I4) C^"^^"(/3)C^"'^°(/4)C"^(|li + I3I). (B.2) 



This statement imphcitly assigns a different role to each leg in the diagram, so each 
permutation of {li, I2, 13, 14} (there are 24), represents in principle a different vertex. 
By analogy to (IB.ip the "4-point vertex-function" of section H] is 



Here (T(li)T(l2)T(l3)T(l4)) includes all the permutations and C{1) are incoming (ACDM 
lensed) propagators. 

B.2 Closing the Diagrams and Symmetry Factors 

The Feynmann rules for the 7 and F vertices include, as usual, division by 2 and 4! 
respectively. This follows from the number of possible permutations of the legs which 
enter the vertex, and is standard in real field theories. The symmetry factor of a diagram 
is the number of different possible ways to generate it. For example since the diagram of 
Fig. [2a] may be generated in two different ways one gets the (1/2)^ x 2 = 1/2 in front of 
the integral (14.61) . 

For the 2-loop calculation we have to draw all possible diagrams which correspond to 
Fig. [2bl For this purpose we have to form loops of the four legs of the F vertex, and then 
add the two 7 vertices in all possible ways. For each diagram of the 24 permutations, 
there are three ways one can form a vaccum energy diagram. These are 



The left diagram is 1-particle reducible, so we ignore it. For each of the other two 
diagrams there are three ways one can insert the two 7 vertices (each of which can be 
inserted in two ways, which cancels the 1/2 that comes with 7). We further note that two 
of the three ways produce topologically identical diagrams, such as 



r(ii, 12,13,14) 



(T(li)T(l2)T(l3)T(l4)) 

cik)c{i2)c{k)c{U) ■ 



(B.3) 
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and therefore all the distinct contributions may be described as 




where we have named the different diagrams, and where each diagram in the top row 
actually represents 2 separate diagrams. 

Therefore, the weights of the diagrams (symmetry factor divided by a factor according 
to the Feynmann rules) are 

48/4! = 2 for each diagram in the top row fB 4) 

24/4! = 1 for each diagram in the bottom row. ^ ' ^ 

B.3 Evaluating the 2-Loop Diagrams 

We now turn to evaluating the 1-particle irreducible 2-loop diagrams. We have 
V^ = w, J ^^^C{h)C{h)C{k)C{k)T{h,-luh^4m-h^ 

©2 = W2 y n^C'(Zi)C'(/2)C'(/3)C(/4)f(li,l2,l3,l4)7(-ll,-l2)7(-l3,-l4)(B.6) 

^3 = j \[-^^C{h)C{h)C{h)C{h)t{\^^^^^ 
f ^ dl 

V, = w, y ]^^C'(Zi)(7(Z2)C'(Z3)C'(Z4)f(li,l2,l3,l4)7(-l2,-l3)7(-li,-l4)(B.8) 
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Figure 5: The four separate non-Gaussian contributions (T'l, V2, T)^ and relative to 
the 1-loop contribution to the S/Nxemp, for the toy model discussed in section HI per /. 
As may be read from the graph, the Pi contribution is positive, while the other three 
contributions are negative. 



respectively, where Wi is the corresponding weight, according to (1B.4I) . 

Applying the toy model ( ]4.ip gives the four separate contributions which may be seen 
in Fig. |5l Simplifying the expressions as much as possible one gets 



V2 



wiA Sip (lo) 
-W2A5ip{\o 



wsASij (lo) 



dl |7(l,lo-l)|^ 
(27r)2 C{1) C(/)C(|lo-l|) 
dl 



dk [l-kfC^ik) 



C"'^i^'^(OC~(|lo + l|) 



unlen / 



(2 



TT 



(27r)2 C(|k 
7(1, -lo - 1) 



II 



(B.9) 



. 1] [k . (lo + 1)] 



C(/)C(|lo + l|) 
1-k) 



(2vr) 



C(|k-l|)C(|lo-li|) 



(B.IO) 



dl C"°''^°(/) |7(l,lo-l)|^ 
(27r)2 C{1) Cil)C{\\o-\\) 

(2vr)2^^ IJi^ C(|k-l|) 
dl 



{2n) 



7(1, -lo - 1) 
C{l)C{l) 



dk 



[k-l] [k-(lo + l 



C"-i-(|lo + l-k|) 



C'^(A;)7(k-l,lo + l-k) 



(B.ll) 

k)] 
(B.12) 



C(|k-l|)C(|lo + l-k| 
The overall contribution of the 2-loop diagrams to the S/N is negative, and yields the 
plateau in the accumulated S/Nxemp as seen in Fig. 
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